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The measurement of peak laser intensities exceeding 10^'' W/cm^ is in general a very challenging task. We 
suggest a simple method to accurately measure such high intensities up to about lO^^W/cm^, by colliding 
a beam of ultrarelativistic electrons with the laser pulse. The method exploits the high directionality of the 
radiation emitted by ultrarelativistic electrons via nonlinear Thomson scattering. Initial electron energies well 
within the reach of laser wake-field accelerators are required, allowing in principle for an all-optical setup. 
Accuracies of the order of 10% are theoretically envisaged. © 2012 Optical Society of America 
OCIS codes: 190.1900,020.2649 
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Since the successful application of the chirped pulse 
amplification technique to generate short laser pulses [T] , 
records in terms of peak intensity are being continu- 
ously set. Intensities as high as 2 x 10^^ W/cm^ have 
already been reached and much higher ones are en- 
visaged [3|, allowing for investigations in different fields 
including strong-field quantum electrodynamics (QED) 
and plasma physics [3] . 

The analysis of experiments employing ultrarelativis- 
tic optical laser pulses, i.e. optical pulses with intensities 
exceeding 10^° W/cm^, requires the precise knowledge 
of quantities characterizing the pulse itself such as its 
peak intensity. However, ultrarelativistic peak intensity 
measurements are especially difficult because the dam- 
age threshold of the equipment is usually far exceeded 
by such strong laser pulses. 

To the best of our knowledge, the peak intensity 
of ultrarelativistic pulses is currently determined by 
measuring energy, duration and spot-size of the pulse, 
the latter however at a lower intensity. Therefore, this 
method is prone to errors since the spot-size can be af- 
fected by increasing the pulse intensity pll5]l6] . It is desir- 
able to find precise methods in order to measure directly 
the peak intensity of the laser field. There are several 
suggestions for direct in situ measurements of peak in- 
tensities by employing multiply-charged ions. The gen- 
eral idea behind these methods is to place ions at the 
laser focus, with charge number selected according to 
the expected intensity. Then, by measuring either the 
photo-ion momentum distribution [7], or the ionization 
fraction [5J|H], one can determine the laser peak inten- 
sity. An alternative method is based on the emission of 
electrons "born" via ionization within the laser pulse [9] . 

In this Letter we suggest a novel method to determine 
the peak intensity of an ultrarelativistic laser pulse, by 
measuring the angular aperture of the radiation spec- 
trum emitted via nonlinear Thomson scattering by an 
ultrarelativistic electron beam colliding head-on with the 
pulse. The method is practically insensitive to the precise 



temporal shape of the laser pulse and it allows for single- 
shot intensity measurements with theoretical accuracies 
in principle of the order of 10%. Units with c = h = 1 and 
a metric with signature H are employed through- 
out. 

Classically it is known that an accelerated ultrarela- 
tivistic electron (mass m and charge e < 0) emits radia- 
tion within a narrow cone around its instantaneous veloc- 
ity, with aperture of the order of l/j{t), with j{t) » 1 
being its instantaneous Lorentz factor [TU]. We exploit 
the high directionality of the radiation emitted by an 
ultrarelativistic electron inside a strong laser pulse to 
relate the angular aperture of the emitted spectrum to 
the peak intensity of the driving pulse. In order to obtain 
an analytical formula we assume that the laser field can 
be approximated as a plane wave. Our numerical results 
will show that this approximation gives accurate results 
also for a tightly focused Gaussian beam. For the sake of 
simplicity we also temporarily neglect effects of radiation 
reaction (RR), i.e. of the action of the radiation emitted 
by the electron on the motion of the electron itself [TU] • 

We consider a plane- wave field linearly polarized along 
the X direction, propagating along the positive z direc- 
tion and we describe it via the vector potential A{tp) = 
— (£'o/wo)/(<p)i, where ip = ujoit — z), luq is the central 
angular frequency of the laser, Eq its electric field am- 
plitude and where f{ip) is a shape function such that 
/(O) = and |/(</j)| < 1. From the exact solution of 
the Lorentz equation in such a plane-wave field [TU] in 
the case of an ultrarelativistic electron with initial four- 
momentum p'^iip = 0) = Pq = eo(l, 0, 0, — /3o) with 
Eq = m7o, /3o = (1 — ^/I'oY''^ and 70 ^ 1, we can write 
the tangent of the angle 6{(p) between the instantaneous 
velocity of the electron and the negative z-axis as 



tan9{ip) 



^0 



7ol-eg/2M/47o' 



2 ' 



(1) 



where ^0 = \£\Eo/muJo. Since the angle 0{(p) practically 
coincides with the instantaneous emission angle in the 
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ultrarelativistic limit, Eq. ([T} shows that the maximal 
emission angle 9m is reached when |/(^)| = 1. By intro- 
ducing the peak laser intensity Iq = Eq/Att, we obtain 
the following simple relation between Iq and 9^' 



Jo [10^" W/cm^] = 0.28 



wo[cV]eo[McV] sin 9,, 
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which is independent of the exact form of f{ip). 

As we have mentioned, we have so far neglected RR ef- 
fects, which, in principle, could become important at the 
intensities considered here [11] . In the realm of classical 
electrodynamics such effects can be included by deter- 
mining the electron trajectory via the so-called Landau- 
Lifshitz equation 
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where is the electron's four- velocity, r its proper time, 
F'^'^{x) the background laser field and a = « 1/137 
the fine-structure constant [HIIO]. The Landau-Lifshitz 
equation has been solved exactly for a plane wave [l2] 
and the solution in the above-described scenario provides 



tan 9{lp) 
where 
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with the prime denoting differentiation with respect to 
its argument. By maximizing the expression of 9{ip) from 
Eq. with respect to (p, one obtains the value of 6'™, 
which includes RR effects. Equations (H])-® show that: 
1) RR effects, being mainly a dissipative effect, will in- 
duce a dependence of 9^ on the pulse-shape function 
/((/)); 2) the expression of tan6'(<p) in Eq. Q reduces to 
that in Eq. ([T]) in the formal limit a 0; 3) RR effects 
become important if (47r/3)Q;(a;o/m)7o^o^ ^ 1' with 
N being the number of laser cycles in the pulse, i.e. if 
eo[MeV]/o[1020 W/cm2]7V/u;o[eV] > 1.5 x 10^. Since RR 
effects act mainly as a damping force, which slows down 
the electron, we expect that their inclusion increases the 
predicted value of 9m at a given Iq . This is confirmed by 
the results in Fig. [T] which displays the angle 9{ip) ne- 
glecting (dash-dotted black curve) and including (solid 
red curve) RR effects for /((/?) = siia^ {(p/ 2 N) sin tp and 
for the following numerical parameters: ujq = 1.55 cV, 
^0 = 100 (/o = 4 X 10^2 W/cm2), 70 = 150 (eq = 
77 MeV) and TV = 10 (pulse duration of 27 fs). As ex- 
pected, at the beginning of the motion the two curves 
coincide, but then larger values of 9{(p) are reached if 
RR effects are included. 
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Fig. 1. (Color online) Comparison of the direction of 
propagation of the electron with (solid red curve) and 
without (dash-dotted black curve) taking into account 
RR, as a function of the plane- wave phase if. Numerical 
parameters are given in the text. 

Below, we present a numerical example showing the 
virtues of our method. We consider a Gaussian beam 
[13] , spatially focused at the origin of the coordinate sys- 
tem, linearly polarized along the x direction, propagating 
along the positive z direction and with a sin^ time enve- 
lope initially centered at —NXo/2, with Aq = 27r/wo be- 
ing the laser central wavelength. The numerical param- 
eters characterizing the beam are: waist size of 5 fim, 
ujQ 1.55 eV, Jo = 2 X 10^^ W/cm^ (^0 = 69) and 
N = 10. Also, the electron is initially at the origin and 
counterpropagates with respect to the laser beam with 
an energy of £0 = 23 MeV (70 = 45). The inner plot in 
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Fig. 2. (Color online) Angular-resolved emitted energy 
spectrum cPW / dfldco (inner plot) and total energy emit- 
ted dW/dil per unit solid angle (outer plot), calculated 
for numerical parameters given in the text. Vertical red 
lines indicate our theoretical predictions for the maximal 
emission angle 9m , with (solid line) and without (dashed 
line) RR included. 

Fig. [2] shows the single-electron emission spectrum [T4] 



d^W 



dtnx{nx /3(i))e*'^(*-"-''(*» 



, (7) 



dfldu 47r2 

i.e. the energy W radiated per unit solid angle and 
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per unit frequency w, on the x-z plane, while the outer 
plot displays the corresponding integrated emitted en- 
ergy dW/dn = dujd^W/dndio. In Eq. Q r{t) and 
f3{t) ~ dr{t)/dt are the instantaneous position and ve- 
locity of the electron calculated via Eq. ([3]) , and n is the 
observation direction. The two vertical red lines indi- 
cate our theoretical expectations, one (dashed) without 
and one (solid) with RR taken into account. The figure 
shows a very good agreement between numerical and an- 
alytical results, although the latter have been calculated 
in the plane-wave approximation. Also, the inclusion of 
RR only increased the predicted value of 9^ by about 
2.2°. We note that the energy required for the electron 
is well within the reach of presently available table-top 
laser-plasma accelerators (LPAs) [TS]. The main advan- 
tages of using a LPA are the microscopic dimensions of 
the electron bunches (transverse radius and longitudinal 
length ~ 10 /^m |15j ) permitting an efficient overlap- 
ping with the laser beam and the allowance for a pre- 
cise synchronization of the electron and the laser beams. 
The typical length of the electron bunches generated via 
LPAs implies that the emission in the relevant part of 
the spectrum (at energies > 10^ eV, i.e. at wavelengths 
smaller than 0.1 nm, see Fig. [2]) occurs incoherently. For 
this reason, in order to simulate an electron beam, we 
have performed simulations by varying the initial con- 
ditions of the electron. We have seen, for example, that 
by shifting the initial transverse electron position within 
a radius of about 1 /^m from the origin, the angular ex- 
tension of the spectrum in Fig. [5] remained practically 
unchanged, while at larger radii it becomes narrower. 
Thus, by including conservatively only those electrons 
in the bunch within a radius of 1 /xm and by assuming 
about 10^ electrons per bunch [15], we can estimate that 
in a single shot about 10^ would contribute to the spec- 
trum as shown in Fig. [2] By estimating the azimuthal 
angle (/> as 1 /70 , we conclude from the numerical results 
in Fig.[2]that about 10^ photons are expected to be emit- 
ted in a single shot at angles between, for example, 73° 
and 75°, which are sufficient for a single-shot measure- 
ment of the intensity. In this estimation we have used 
~ 10^ eV as average emitted photon energy. 

By indicating as /S.F the uncertainty on the generic 
quantity F, Eq. ([2]) implies that 




showing that at 70 w ^0 /2 (so that sin ss 1) the contri- 
bution of the uncertainty in 6m is minimal. By safely ne- 
glecting Awq / ljJq and by assuming Aeq / £0 ~ 5% }15j and 
A^m « 1° (corresponding to the ratio between the di- 
mension of a typical high-energy photon detector 1 cm 
and the distance ~ 1 m of the detector from the inter- 
action region), we obtain AIq/Iq ~ 10% in our numeri- 
cal example. This theoretical estimate would in principle 
significantly improve the accuracy of intensity measure- 
ments, which nowadays are often even larger than 50%. 



In the above analysis we have neglected quantum 
effects, which mainly amount in the ultrarelativistic 
regime to the photon recoil experienced by the electron 
in the photon emission [14]. These effects become im- 
portant at X = 27o-E/-Ecr ^ 1, with Ecr = iv?/\e\ = 

I. 3 X 10^6 V/cm being the "critical" field of QED [l^ 
and in our numerical example it was x ^ 2 x 10~^. 
In the convenient regime 70 ^ Co/2 quantum effects 
can be neglected {x < 0.1) for intensities smaller than 

10^3 W/cm^ at Wo = 1.55 eV, which determines the 
intensity upper limit of our method at least for an op- 
timal accuracy. In fact, at optical intensities larger than 
10^2 W/cm^ an electron emits on average more than one 
photon already in one laser period [T4| . The classical ex- 
pression of the emitted spectrum in Eq. ([7]), which is 
valid at X ^ 1. automatically takes into account multi- 
ple photon emission |16j . However, at x ^ li when 
recoil effects become important, there are not yet QED 
calculations on angular-resolved multiple-photon spectra 
in a strong laser field. 
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